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ABSTRACT 

We present a sample of radio afterglow light curves with measured decay 
slopes which show evidence for a flattening at late times compared to optical 
and X-ray decay indices. The simplest origin for this behavior is that the change 
in slope is due to a jet-like outflow making a transition to sub-relativistic ex- 
pansion. This can explain the late-time radio light curves for many but not all 
of the bursts in the sample. We investigate several possible modifications to 
the standard fireball model which can flatten late-time light curves. Changes 
to the shock microphysics which govern particle acceleration, or energy injection 
to the shock (either radially or azimuthally) can reproduce the observed behav- 
ior. Distinguishing between these different possibilities will require simultaneous 
optical/radio monitoring of afterglows at late times. 

Subject headings: gamma rays:bursts — radio continuum: general — ISM:jets and 
outflows — shock waves 



1. Introduction 

One of the defining characteristics of X-ray and optical afterglows is the observed power- 
law decay of their light curves. On timescales of hours to days after a burst the exponent 
a (defined by F^, oc typically lies in the range of —1 to —2. The light curves of several 
GRBs have also been seen to undergo an achromatic break, steepening their temporal decay 
by Aa ~ 1. The standard afterglow model provides a framework in which to interpret 
these power-law slopes and their changes. Synchrotron emission is produced in a relativistic 
shock which accelerates electrons to a power-law distribution with energy index p given by 
N(7e) oc above some minimum energy 7^. The evolution of the expanding blast wave 
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is sensitive to the geometry of the shock, the kinetic energy in the shock, and the structure 
of the circumburst medium. This leads to the temporal (and spectral) evolution of the 
light curves whose power-law indices are predicted for different cases (Sari, Piran & Narayan 
1998; Sari, Piran & Halpern 1999; Chevalier & Li 2000). The most commonly accepted 
explanation for the sharp breaks in optical and X-ray light curves is that the outflows are 
coUimated and that the change in a is the result of the Lorentz factor F dropping below 
the inverse opening angle of the jet (Rhoads 1999). However, it should be noted that 
the jet signature is not identified unambiguously in all bursts, and for some events there are 
alternate explanations for the origin of the break (Kumar & Panaitescu 2000; Masetti et al. 
2001; Wei & Lu 2002). 

The situation at radio wavelengths is slightly more complex. Initially, at least, the 
emission comes from the lowest energy electrons at which are radiating only a small 
fraction of their energy at radio frequencies below i/^- This gives rise to the familiar spectral 
slopes of u^^^ and in the optically thin and thick part of the radio spectrum, respectively, 
and it produces light curves that rise with time. A jet break is expected initially to produce 
only a shallow power-law decay (e.g. t^^^'^ to t") of the radio light curve until a time when 
the synchrotron peak z/^ passes through the band. If the expansion remains relativistic, the 
post-jet power-law decay aji will be the same as that at optical and X-ray wavelengths {i.e., 
a ~ -2) 

In addition to the late onset of the steep jet decay, it may be possible to use the long- 
lived nature of the radio emission to look for changes in a when the expanding shock slows to 
subrelativistic speeds (Waxman, Kulkarni & Frail 1998). Deviations in the power-law decay 
could also be produced by density enhancements in the circumburst medium on parsec scales, 
similar to that seen in supernovae (Montes et al. 2000). Likewise, it is possible that the 
bulk of the GRB blast energy is in a low F component and that the late-time afterglow 
will be refreshed as this slower shock catches up to the deccelerating swept-up shell. There 
may also exist an additional low energy electron component (Waxman 1997; Bhattacharya 
2001; Li & Chevalier 2001; Panaitescu 2001), which while energetically unimportant it could 
significantly modify the late-time afterglow light curves. Finally, all of these processes could 
be masked by the emergence of an underlying host galaxy, whose radio emission is the result 
of prodigious star formation. 

Motivated by the possibility of detecting these effects at late times, we have analyzed a 
set of well-sampled GRB radio light curves for which it is possible to measure their temporal 
decay indices. In §2 we compare the sample of X-ray and optical afterglows, while in §3 
the comparison is carried out with the radio and optical afterglow sample. In §4 we discuss 
known sources of biases that affect decay slope measurements, while in §5 and §6 we discuss 
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possible explanations for the observed flattening of the radio hght curves. 

2. An X-ray/Optical Compcirison 

In Table 1 we provide a list (complete until 2002 January) of all afterglows with ac- 
curately measured temporal decay indices at both X-ray and optical wavelengths. Most of 
the X-ray measurements were obtained between 8 and 24 hours after the burst. For the 
optical data we list the time interval over which the power-law flt was made. If a jet break 
occurred in the light curve, two values of cto are given along with the relevant timescale. 
When required, a constant emission component from a host galaxy was also fit to the optical 
light curves. 

In the top panel of Figure 1 we show a plot of versus cKq. Most of the points are 
distributed about the line = o^o with small scatter but there are a few significant deviations 
where the X-ray decay is steeper than the optical decay. For reasonable physical parameters 
for the shock in a constant density medium, the synchrotron cooling frequency Uc is expected 
to lie between the optical and X-ray bands on these timescales (Panaitescu & Kumar 2001). 
The temporal slope of light curves measured above z/g will be steeper by an amount Aq;=1/4 
(Sari et al. 1998), consistent with all points which he above the fudicial ax = olo fine. 

More than half of the points in Fig. 1 have a tendency to lie below the ax = ao line. In 
the simplest afterglow models ao is not expected to be steeper than ax- We could expect such 
a result if the shock wave were propagating into a density gradient, such as that produced 
by mass loss from a progenitor star. For reasonable physical parameters Vc is expected to 
lie below the optical band for the timescales of interest here (Li & Chevaher 2001). A more 
plausible explanation is that inverse Compton emission, which has been reported for several 
bursts (e.g., Harrison et al. 2001), is producing a slight flattening the X-ray light curves. 
Finally, for GRB 980519 it is likely that the jet break occurred close to the time of the X-ray 
measurements (Nicastro et al. 1999; Jaunscn et al. 2001) producing an estimate of ax that 
is intermediate between the two asymptotic values. 

3. A Radio/Optical Comparison 

The 10 afterglows in Table 1 were selected from Frail et al. (2003) to have well-sampled 
light curves at a frequency of 8.5 GHz. The remaining 15 GRBs in this catalog do not have 
sufficient measurements to determine a decay index aji at late times. A least-squares fit was 
made to each dataset. In order to be as free as possible from interstellar scintillation and 
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deviations from power-law decay due to curvature effects in tfie syncfirotron spectrum (see 
§4), each burst was fit to several different starting epochs. The final values of listed in 
Table 1 were often a compromise between making the starting epoch as late as possible while 
including sufficient data points for a fit. In Figure 2 we show the results of our fits for the 
four best sampled events. The reduced ranges from 0.6 to 1.3, typical of the sample as a 
whole. 

In the bottom panel of Figure 1 we show a plot of versus Oq. The optical decay 
index was chosen over for this comparison because there are more joint radio/optical 
datascts and because their lightcurvcs arc taken over similar timescales. Nonetheless, for 
the discussion that follows it is important to realize that the overlap between the optical 
and radio time intervals is relatively small. This will complicate attempts to understand the 
relation between and (see §6). 

There is a clear trend in Figure 1, namely the radio decay indices are equal to or 
substantially flatter than the optical (or X-ray) decay indices. Moreover, the an values 
rarely exceed —1 and there are no examples of ~ —2. The data can be grouped into 
three different regions in Figure 1. The first are those with an ~ cto, the second have 
CKo < —2, and the third group have cto > —2. In the next several sections we discuss several 
possible explanations for the origin of these temporal slopes. 



4. Systematic Sources of Light Curve Flattening 

Some of the measured ag values are not as steep as is expected for a jet break due to the 
presence of a host galaxy, which for most events dominates the optical light curves typically 
between one week and one month after the burst. In those cases, contamination from the 
host galaxy may prevent a true determination of the post-jet Radio data have been used 
in several cases (e.g., GRB 000418 and GRB 980703) to further substantiate claims of jet 
breaks (Berger et al. 2001; Frail et al. 2003) but removing the host galaxy contribution is 
often difficult. A second source of fiattening for optical light curves may be the excess of 
optical fiux at ~20(l+2;) days, commonly attributed to the rise of a supernova component 
(e.g. Bloom et al. 2002). 

A host galaxy may also fiatten the late-time radio light curves. For GRB 980703 
there is good evidence that the fit requires an additional contribution from an underly- 
ing galaxy (Berger, Kulkarni & Frail 2001). A starburst host galaxy has also been proposed 
for GRB 000418 (Berger et al. 2003). However, while the submillimeter detection of the 
host is robust, in our view further multi-frequency radio observations are required before the 
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centimeter identification of a host can be considered certain. Adding a host component to 
the GRB 000418 fit yields aR = -1.72 ± 0.09 and fhost^^7 A*Jy- 

We do not expect radio galaxies to be detected as frequently as those at optical wave- 
lengths. In order to significantly flatten the values in Table 1 we require Fhost ~ 30 — 40 
//Jy, which at 8.5 GHz for the typical redshifts of these events, requires SRF~ 10^ Mq yr~^, 
or Lboi > 10^^ L0. If GRBs trace star formation in the universe, as is currently believed, 
then only about 20% of GRBs are expected to be found in such ultraluminous host galaxies 
(Barnard et al. 2003). This fraction is in good agreement with a more direct determi- 
nation from centimeter and submillimeter observations of GRB host galaxies (Berger et al. 
2003). Thus we expect that one, perhaps two, afterglows in our sample have undetected host 
galaxies. Fortunately, continued multi-frequency radio monitoring can be used to detect and 
subtract off any galaxy contribution. 

A second, more insidious, problem is to ensure that ur is not measured when the light 
curve is in transition between two power-law behaviors. For example, in the case a jet-like 
outflow there is a timescale tj when the edge of the jet becomes visible and the optical 
emission steepens to cKq = —p- In contrast at this time, the radio flux (which is typically 
emitted below the synchrotron peak Vm at this time) decays with ur — —1/3 until Vm 
passes through the radio band, after which it decays as ur = —p in the relativistic phase 
or as aR = (21 — 15p)/10 in the sub-relativistic phase. Fits which encompass any of these 
transitions will give values of aR which are in between these extremes. 

An artificial flattening of the light curve will be measured if the fit is made before i/m 

passes through the radio band. If Um passes through the optical band at a time to it is 
straightforward to show that for a jet Vm passes through 8.5 GHz at tR ~ 230t]/^ x t^^/^ 
days provided to < tj < tR. Thus for typical shock parameters we expect t/j ~ 5 — 30 days 
(Sari et al. 1999). Most of the fits in Table 1 were made at a time t > tR but we caution the 
reader that the shallow aR values for at least two GRBs may be the result of fitting over 
this transition. The only way to improve on these fits is to reduce the dependence on earlier 
measurements through deep observations of the faint afterglow emission at late times. 

5. A Dynamical Origin for the Late Decay Slopes 

A change in the temporal slope is expected in the basic afterglow model when the 
expansion of the blast wave becomes subrelativistic. This dynamical transition has been 
predicted for some time (Wijers, Rees & Meszaros 1997), and it has been claimed to have 
been seen in a number of events (Frail, Waxman & Kulkarni 2000; Dai & Lu 2000; Piro et 
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al. 2001, Int'Zand et al. 2001). This occurs on a timescale when the rest mass energy swept 
up by the expanding shock becomes comparable to the initial kinetic energy of the ejecta. 
For kinetic energies of 10^^ and circumburst densities of 1 cm~^ this occurs on a timescale 
of order 100 days. After this time the dynamical evolution of the shock is described by 
the Sedov- Taylor solutions rather than the relativistic formulation of Blandford & McKee 
(1976). 

Independent of geometry, the expected temporal slope in the non-relativistic regime is 
(^NR = (21 — 15p)/10 for a constant density (ISM) medium, and a^R = {5 — 7p)/6 for a wind- 
blown medium {i.e., p oc r~^). These values assume that the synchrotron break frequency 
Um has passed through the band but the cooling frequency Uc remains above the band. Livio 
& Waxman (2000) provide a convenient table of a's for different cases. 

For a spherical fireball undergoing a transition to non-relativistic expansion, the light 
curves are expected to steepenhy Act ~ 0.3, while for a jet-like expansion the light curves are 
expected to flatten by Aa ~ 1. By and large the aR values in Table 1 are flatter than the ao 
(and ax) values measured at earlier times. This immediately rules out a spherical geometry 
for the majority of bursts in Table 1 with measured aR values. For a jet-like outflow it is 
expected that ao — —p (for t > tjet) and therefore the values of aNR can be calculated and 
compared to aR. For bursts with ~ -2.2 (GRB 991208, GRB00031C, GRB 000418 and 
GRB 000926) we expect a^R ~ —1.2 in the ISM model and a^R ~ —1.7 in the wind model. 
From Fig. 1 we note that at least for these GRBs the JET+ISM model provides a better 
(but not ideal) fit to the observed aR values than a JET+WIND model. The deviations are 
in the sense that the observed aR are systematically more shallow than expected from the ao 
measurements. This is not likely to be a systematic bias because it is in the opposite sense 
of what would be expected if the optical decay was still in transition to its asymtoptic value 
{i.e., ao > -p). The JET-MSM model may also explain events like GRB 980703 for which 
there is good evidence that the temporal slope ao is underestimated because the optical light 
curves are contaminated by an underlying host galaxy (§4). 

A simple jet model does not work for GRB 970508, for which ctj. > «(, ~ aR. We note 
that both the optical and radio temporal slopes were measured at comparatively late times 
and therefore both of the fits may be dominated by points when the afterglow was in the 
sub-relativistic phase. In this case a near-spherical shock with p ~ 2.3 provides a good 
description of the light curves (Frail et al. 2000, but see Chevalier & Li 2000). The most 
difficult challenge for the JET+ISM model comes from the third cluster of GRBs in Fig. 1 
with ao > —2. In two cases (GRB 010222, and GRB 000911) adopting tto — P yields estimates 
of aNR that are substantially flatter than the observed aR values. The JET-I-WIND model 
works better but it is not clear that this model can explain the pre-jet behavior of the light 
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curves. 

To summarize, while a sub-relativistic transition of a jet-like outflow in a constant 
density medium provides a reasonable explanation for the flattening of the radio light curves 
for some of the GRBs signiflcant departures from this behavior are seen (e.g. GRB 000926). 
The radio slopes are substantially flatter than expected in this simple model and some other 
source must be found. 



6. Alternate Origins and Conclusions 

In the previous sections we have presented evidence that the temporal decay of radio 
light curves measured months after the burst show evidence for a flattening compared to 
the decay slopes measured at optical and X-ray wavelengths within the flrst week after the 
burst. The simplest explanation, consistent with most of the data, is that the flattening is 
the result of a dynamical transition of the shock to sub-relativistic expansion. 

There are other physical effects that may lead to the flattening of radio hghtcurves. 
Several authors (Bhattacharya 2001; Li & Chevalier 2000; Panaitcscu 2001; Dado, Dar & De 
Riijula 2003) have modified the electron energy distribution by introducing a break in the 
power-law below which the spectrum is hard {i.e., p < 2). This has the virtue of requiring 
no additional source of energy while simultaneously fitting for the different decays slopes in 
the X-ray, optical and radio bands. Its disadvantages are that it requires the introduction of 
another free parameter in the modeling, and that current simulations of particle acceleration 
in ultrarelativistic shocks are unable to produce hard energy spectra (Achterberg et al. 
2001). Other modifications to the shock microphysics are possible (Rossi & Rees 2003). For 
example, Yost et al. (2003) have relaxed the usual assumption that the magnetic energy 
density behind the shock is constant, and note that the late-time light curves flatten when 
the magnetic energy grows inversely with the Lorentz factor of the shock. 

The flattening could also be maintained by a continuous or episodic injection of energy, 
rather than the one-time injection as is commonly assumed (Rees & Meszaros 1998; Sari & 
Meszaros 2000; Kumar & Piran 2000). Slower moving shells of ejecta catch up to the dec- 
celerating main shock and re-energize it, causing to afterglow to brighten at all wavelengths 
(e.g. Panaitescu, Meszaros & Rees 1998). This "refreshed shock" model has been invoked to 
explain the optical behavior of GRB 010222 (Bjornsson et al. 2002) and GRB 021004 (Fox 
et al. 2003). This explanation is problematic for the radio flattening because it requires 
continuous energy injection but with a delayed turn-on in order to both maintain the shal- 
low decay of the radio light curve, while preserving the steeper optical decay. One way to 
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overcome this difficulty would be to add energy though a two component jet-like outflow, 
with the radio emission originating from an outflow with a wide opening angle carrying the 
bulk of the energy. This was flrst proposed for GRB 991216 (Frail et al. 2000), but Berger et 
al. (2003) recently have made a stronger case for GRB 030329. Since the bulk of the energy 
in this case is carried by the slower moving cjccta a crucial test would be to carry out a 
minimum energy analysis like for GRB 970508 (Frail et al. 2000) and look for large excesses 
compared to energies derived from conventional methods (Panaitescu & Kumar 2002; Frail 
et al. 2001). 

In principle it should be possible to distinguish between these alternate explanations 
for the observed flattening. A dynamical transition is achromatic, and so a break in the 
optical and radio light curves should occur at the same time. A modified electron spectrum 
will produce two spectral components, recognized by comparing the spectral slopes both 
within and between the radio and optical bands. Energy injection, added either radially via 
refreshed shocks, or azimuthally from complex jet structure, will most likely be identified 
from multi-component light curves. All of these tests require near-simultaneous optical and 
radio light curves but as is evident from Table 1 such data is currently lacking. Although 
there are several practical problems to overcome (§4), future monitoring of optical afterglows 
should be extended to produce a better overlap with the late-time radio measurements. The 
recent bright and nearby GRB 030329 is a promising candidate. 

The National Radio Astronomy Observatory is a facility of the National Science Foun- 
dation operated under cooperative agreement by Associated Universities, Inc. DAF thanks 
the Astronomical Institute at the University of Amsterdam for their hospitality during the 
time when this paper was being written. 
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Table 1. X-ray, Optical and Radio Decay Slopes 









Epoch 




Epoch 


Ref. 


GRB 




Oio 


(days) 




(days) 




970228 


-1.33 ±0.12 


-1.58 ±0.28 


1-5 






1,2 


970508 


-1.1 ±0.1" 


-1.30 ±0.05 


2-120 


-1.34 ±0.10 


115-309 


3,4 


971214 


-0.96 ±0.15 


-1.20 ±0.02 


0.5-3 






5,6 


980329 


-1.35 ±0.03 


-1.21 ±0.13 


0.7-10 


-1.15 ±0.17 


55-121 


7,8 


980519 


-1.83 ±0.3 


-2.05 ± 0.04 


0.35-2 






9,10 


980703 


-1.24 ±0.18 


-1.61 ±0.12 


0.8-10 


-1.33 ±0.06^ 


27-210 


5,11 


9901/0 


— 1.41 ± 0.05 


— 1.10 ± 0.0.5 


0.18-2 






5,12 


990510 


-1.0 ±0.1 


-0.82 ±0.02 


0.15-1.2 






13,14 


991208 




-2.2 ±0.1 


2-7 


-1.07 ±0.09 


53-293 


15 


991216 


-1.61 ±0.06 


-1.22 ±0.04 


0.4-2 






16,17 






-1.80 ±0.30 


2-15 


-0.85 ±0.16 


8-78 




000301C 




-2.29 ±0.17 


4-13 


-0.93 ±0.12 


45-165 


18 


000418 




-1.41 ±0.08'^ 


2-14 


-1.05 ±0.10 


75-202 


19 


000911 




-1.46 ±0.05 


1-15 


-0.91 ±0.13 


3-23 


20 


000926 


-1.89 ±0.18 


-2.38 ±0.07 


2-30 


-0.76 ± 0.08 


25-288 


21,22 


010222 


-1.33 ±0.04 


-1.57 ±0.04 


1-50 


-0.55 ± 0.09 


5-206 


23,24 



Note. — The columns are (left to right): (1) GRB name, (2) the X-ray temporal decay 
index defined by Fx oc f^^ , (3) the optical temporal decay index, (4) the timerange over 
which the power-law fit was made to the optical light curves, (5) the radio temporal 
decay index, (6) the fit timerange, (7) references for the X-ray and optical decay slopes. 
" The X-ray decay index for GRB 970508 is approximate since there was substantial 
temporal and spectral variability for this burst. ^ Based on the findings of Berger, Frail 
& Kulkarni (2001), we subtracted a constant component for the host galaxy of this 
burst. Berger et al. (2001) find evidence for a jet break at t~26 d based on radio 
and late-time optical measurements. A fit to the data with a standard jet model in a 
constant density medium gives reasonable solutions for ag — 2.4. 

References. — (1) Costa et al. (1997); (2) Reichart (1999); (3) Pruchter et al. (2000); 
(4) Piro et al. (1999); (5) Stratta et al. (2001); (6) Diercks et al. (1998); (7) in 'T 
Zand et al. (1998); (8) Reichart et al. (1999); (9) Nicastro et al. (1999); (10) Halpern 
et al. (1999); (11) Vreeswijk et al. (1999); (12) Kulkarni et al. (1999); (13) Plan et al. 
(2001); (14) Harrison et al. (1999); (15) Sagar et al. (2000); (16) Frail et al. (2000); 
(17) Halpern et al. (2000); (18) Jensen et al. (2001); (19) Berger et al. (2001); (20) 
Price et al. (2002); (21) Piro et al. (2001); (22) Harrison et al. (2001); (23) in' t Zand 
et al. (2001); (24) Galama et al. (2003). 
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Fig. 1. — The top panel shows the decay indices a defined by F,^ oc for a sample of 
gamma-ray bursts with X-ray and optical afterglows. The solid line corresponds to ax = ag, 
while the dashed lines have constant offsets of ±0.25. The bottom panel shows decay indices 
for a sample of bursts with radio and optical afterglows. The solid line is — ag while the 
dotted and dashed lines are the expected relation between the relativistic and non-relativistic 
temporal decay slopes for a jet-like outflow in constant density and wind circumburst media, 
respectively. In plotting these lines it is assumed that cKq = —p and a.NR — ocr. 
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Fig. 2. — Radio light curves of four GRBs at a frequency of 8.46 GHz. A least squares fit 
was made to the decaying portion of the light curve. The time interval and the slope of 
each fit is indicated by the solid lines. For GRB 980703 a constant flux density was added to 
the fit to account for the emission from an underlying host galaxy. For GRB 000418 a pure 
power-law fit is shown along with a host component included. The slopes of the optical light 
curves and the time interval over which the fits were made are illustrated schematically with 
dashed hues. 



